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Using the density functional theory (DFT) with the generalized gradient approximation (GGA), 
the structural and electronic properties of wurtzite A1N, GaN, InN, and their related alloys, 
ALGai-xN and In^Gai-^N, were calculated. We have performed accurate ab initio total energy 
calculations using the full-potential linearized augmented plane wave (FP-LAPW) method to in- 
vestigate the structural and electronic properties. In both alloys we found that the fundamental 
parameters do not follow Vegard's law. The lattice parameters, a, c, and u, for the AUGai-^N 
alloy are found to exhibit downward bowing, while for In^Gai-^N there is an upward bowing for 
the a and c parameters and a downward bowing for the internal parameter, u. Furthermore, we 
found that for both alloys, the band gap value does not follow Vegard's law. As a by-product of our 
electronic band structure calculations, the effective masses of the binary compounds as well as their 
related alloys were calculated. We show that the calculated properties for the binary compounds, 
as well as for the studied alloys, show good agreement with most of the previously reported results. 
Finally, using the frozen phonon approach, the Ai (TO) mode for the different systems studied in 
this work was calculated. Our calculations show good agreement with experimental values reported 
for the binary compounds. For the ternary alloys, our calculations reproduce experimental values 
for Al^Gai-^N as well as theoretical predictions for In^Gai-^N. 
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I. INTRODUCTION 

The Ill-Nitride semiconductors have attracted much 
attention over recent years because of their potential ap- 
plications in technological devices. This is due mainly 
to the fact that the energy gap can be tuned over a 
wide spectral range from the visible to the ultraviolet 
regime of the electromagnetic spectrum. Although the 
zinc-blende and wurtzite structures are present in the 
GaN, A1N, and InN semiconductors, it has been demon- 
strated experimentally that wurtzite is the most stable 
structural phase of these compounds. Moreover, due to 
their high chemical and thermal stability, the III-Nitridcs 
are ideal candidates for applications under extreme con- 
ditions such as high temperature applications. In the 
wurtzite crystalline structure the value of the band gap 
ranges from 0.8 eV for InN 3.4 eV for GaN Q and 6.2 
eV for A1N [3j , providing a huge interval of energies for 
this parameter whenever the concentration forming the 
alloy is carefully selected. The hexagonal wurtzite struc- 
ture is extensively utilized because all of the Ill-nitride 
semiconductors and their related alloys exhibit a direct 
band gap energy, which results in a high emitting perfor- 
mance [j, [f| . Due to the remarkable progress in epitaxial 
growth technology, high quality samples of these com- 
pounds can be produced. High-quality wurtzite InN is 
currently available and its direct band gap energy has 
been determined to be around 0.8 eV, which is much 
smaller than the commonly accepted value of 1.9 eV 
From a theoretical point of view, many calculations us- 
ing different methods have been done to characterize the 



structural, electronic, and optical properties of these sys- 
tems; however there is still no agreement in the scientific 
community concerning the values of certain parameters, 
since they show significant scattering when we compare 
the experimental or theoretical results published in liter- 
ature, as we discuss in this work. 

In this paper, by means of numerical calculations based 
on first principles, we present a study of the struc- 
tural and electronic properties of the A1N, GaN, InN 
semiconductors and their related alloys, A^Gai-aN and 
In a; Gai_ a ;N, with a wurtzite structure. The analysis was 
made by calculating the total energy from first principles. 
First, we analyzed the binary compounds, GaN, A1N, 
and InN, and then their related alloys, Al^Gai-^N and 
In^Gai-^N, for x = 0, 0.25, 0.50, 0.75, and 1.0. Our 
calculations were based on the density functional the- 
ory (DFT) using the generalized gradient approximation 
(GGA) in order to calculate the exchange-correlation 
term in the total energy. We used the Wien2k simu- 
lation package developed by the Vienna University of 
Technology. In section II, we describe the model used 
in this work, while in section III, we discuss our results 
and compare them with data found in related literature. 
Finally, in section IV, we present our conclusions. 



II. THEORETICAL ASPECTS AND 
COMPUTATIONAL METHOD 

Our calculations were performed within the framework 
of the density functional theory (DFT) Q , which states 
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that all of the ground state properties of a system are 
functionals of the electron density and the total energy 
is expressed in terms of the electron density rather than 
the wave function. At present, DFT is one of the most ac- 
curate methods to calculate the structural and electronic 
properties of solids. We have used the full-potential lin- 
earized augmented plane wave method (FP-LAPW) as 
implemented in the Wien2k code [3]. As most of the 
first principles methods, LAPW is a procedure used to 
solve the Khon-Sham set of equations for the density of 
the ground state, the total energy, and the eigenvalues 
of a many-electron system. In the present analysis, the 
exchange-correlation energy of the electrons was treated 
using the generalized gradient approximation (GGA) [1]. 
To minimize the energy, the Wien2k code divides the unit 
cell into non-overlapping spheres centered at atomic sites 
(muffin-tin (MT) spheres), and the interstitial region. In 
the MT spheres, the Khon-Sham orbitals are expanded 
as a linear product of radial functions and spherical har- 
monics, and as a plane wave expansion in the interstitial 
region. The basis set inside each MT sphere is divided 
into core and valence subsets. The core states are treated 
within the spherical part of the potential only and are 
assumed to have a spherically symmetric charge density 
that is completely confined within the MT spheres. In 
this work, the valence part was treated within a poten- 
tial that was expanded into harmonics up to I = 10. We 
have used MT sphere radii of 1.6 Bohr for N, 1.9 Bohr 
for Al, 2.0 Bohr for Ga, and 2.33 Bohr for In. The self- 
consistent calculation was considered to converge when 
the total energy of the system was stable to within 10~ 5 
Ry. Care was taken to assure the convergence of the 
total energy in terms of the variational cutoff energy pa- 
rameter. Furthermore, we have used an appropriate set 
of k-points to compute the total energy. To calculate 
the convergence of the total energy we wrote the basis 
functions up to a cutoff radius of R m tK max — 7 Ry for 
both of the binary compounds and the Al^Gai-^N al- 
loy, and R.mtK max = 8 Ry for the In^Gai-^N alloy. We 
have minimized the total energy using different sets of 
k-points in the irreducible part of the Brillouin zone and 
constructing an appropriate grid in the unit cell accord- 
ing to the Monkhorst-Pack procedure @ . The number of 
k-points used was chosen in order to assure convergence 
within our accuracy criterion (10 -5 Ry). 

In the wurtzite structure, the positions of the atoms 
inside the unit cell are (0, 0, 0) and (2/3, 1/3, 1/2) for the 
cation Al, Ga, In; and (0,0, u) and (2/3, 1/3, 1/2 + u) 
for the anion N, where u is the internal parameter for 
the cation-anion separation. We began our study by op- 
timizing the structural parameters for the binary com- 
pounds, GaN, A1N, and InN, starting from the ideal 
wurtzite structure with a ratio c/a = 1.633 and u = 0.375 
for the internal parameter. This optimization was made 
by an iterative process as a function of the volume, V, 
the c/a ratio, and the internal parameter, u, until the 
total energy converged to within 0.01 mRy To model 
the Al x Gai_ x N and In^Gai-^N ternary alloys, we used 



a 32-atom supercell with periodic boundary conditions. 
This corresponds to a 2 x 2 x 2 supercell which is twice 
the size of the primitive wurtzite unit cell in both direc- 
tions: along the basal plane and along the c— axis. We 
minimized the total energy for different values of the con- 
centration, x (0.25, 0.50, and 0.75), as a function of the 
three variables mentioned above. The atomic electronic 
configuration used in our calculations was: Al (Ne, 3p, 
3s), Ga (Ar, 3d, 4s, 4p), In (Kr, 4d, 5s, 5p), N (He, 2s, 
2p). The Ga3d and In4d electrons are treated as va- 
lence band states using the local orbital extension of the 
LAPW method Q. 



III. RESULTS AND DISCUSSION 

A. Structural parameters for the Al^Gai-^N and 
In^Gai-^N alloys 

1. Al x Ga-i- x N 

Table 1 summarizes our results and compares them 
with experimental and theoretical results reported pre- 
viously using different methods. In going from GaN to 
A1N, when the Al-content increases, the values of the lat- 
tice parameters of the ALjGai-^N alloy decrease. This 
is due to the fact that the size of the Al atom is smaller 
than the Ga atom. This is not the case for the inter- 
nal parameter, u, in which we observe an increase of this 
parameter when we increase the Al content of the alloy. 
Figure [1^ depicts the behavior of the a and c parameters 
as a function of the Aluminum concentration. We can see 
from the figure that these parameters show a clear devi- 
ation from the linear behavior stated by Vegard's law, 
which determines the parameter behavior of these alloys 
with a zinc-blende structure [l(J. The deviation from 
Vegard's law can be quantified by adjusting the curves 
in Figure []Jl to the following formula: 

A(x) = xAain + (1 - a;)^GaN - e A x{l - x); (1) 

where A(x) stands for the different structural parame- 
ters, a, c, and u, of the Al x Gai_ x N alloy. Aain (^Gaisr) 
represents the structural parameters of the binary A1N 
(GaN) compound and e A is the respective bowing pa- 
rameter for the lattice and internal parameters. If we fit 
our calculated values from Table U to Eq. ([1]) , we obtain 
a set of bowing parameters: e a = 0.016 A, e c = 0.119 
A, and e u = 0.002. We observe that all bowing param- 
eters have a positive value, which indicates a downward 
bowing, being the lattice constant c which possesses a 
larger deviation from the linear Vegard's law, as was ex- 
perimentall y r eported by S. Yoshida et al. [la ], and by 
Yun et al. [l4j. Other works report that the lattice pa- 
rameters follow Vegard's law (see for example Angerer et 
al. [IB]); however, from from a theoretical approach, this 
alloy has only been studied using the virtual crystal ap- 
proximation (VCA) by M. Goano et al. [ll| or by using 
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first principles by Z. Dridi et al. [l2j]. Other authors have 
reported a non linear behavior of the lattice parameters, 
but their work demonstrates an upward bowing of the 
parameters [lq . 

To compare our results with the available experimen- 
tal information, we plot in Fig. []Jl the calculated value 
of the c parameter obtained by Eq. |T]) together with the 
experimental results found in literature [131 ]. As we can 
see, our calculated values show the same trend depicted 
by the experimental data. Other examples of experimen- 
tal studies for this alloy can be found in the works by D. 
K. Wickenden et al. [l7j], and by K. Itoh et al. pi ]. 



2. In x Ga-i- x N 

The recent developments in blue-green optoelectron- 
ics are essentially due to the high efficiency luminescence 
of In^Gai-zN/GaN heterostructures. Despite their im- 
portance, several properties of In K Gai_ K N alloys are not 
fully understood. For example, the optical properties 
of InN crystals are poorly known since the available 
growth techniques have not allowed the production of 
high quality epitaxial layers. Recent improvements in 
the molecular beam epitaxy (MBE) technique have led 
to the availability of high quality InN films. Photolumi- 
nescence measurements of these films indicate an energy 
gap around 1 eV or less The In x Gai_a;N alloy has 
been studied theoretically by several groups using differ- 
ent methods. M. Goano et al. [TTf] used pseudopotentials 
to compute the gap through the virtual crystal approx- 
imation approach. Z. Dridi et al. [12] used LDA FP 
LAPW and the virtual crystal approximation approach. 
J. Serrano et al. [2(| worked within the framework of the 
density functional theory (DFT) with the local density 
approximation (LDA) using the Ceperley-Alder form for 
the exchange- correlation energy. C. Stampfl et al. [27| 
utilized the DFT, the LDA, and the GGA of Perdew et 
al. for the exchange-correlation functional. A. Zoroddu 
et al. [28| from first principles within the DFT uti- 
lized the plane-wave ultrasoft pseudopotential method 
within both the LDA and the GGA. And P. Carrier et al. 
p^ ] used plane-wave pseudopotentials and the LAPW 
method with the LDA. 

In Table [IT] wc summarize our results and compare 
them with some of the representative theoretical and ex- 
perimental results found in literature. We can see from 
these results that the values of the a and c parameters 
increase when the In concentration increases. We plot 
these results in Fig. [TJd. It is clear from the figure that, 
as in the previous case, there is not a linear dependence 
in these two parameters with an increase in the In con- 
centration. If we adjust these results using Eq. (Q]), as 
we did for the Al^Gai-^N alloy, we obtain: e a = —0.140 
A, e c = -0.188 A, and e u = 0.0001. For this alloy we can 
observe that the deviation parameters for both a and c 
have a negative value, which implies an upward bowing 
and is clearly observed in Fig. QJ>. This is not the case 



for the internal parameter, u, which has a nearly linear 
dependence with the In concentration as demonstrated 
by the bowing parameter, e u = 0.0001. The simulation 
results indicate that the c lattice constant has a larger 
deviation from the linear Vegard's law when compared 
with the lattice constant a. Finally in Fig. [21 we plot the 
internal parameter, u, for both alloys as a function of the 
concentration, x. Solid circles and squares correspond to 
our theoretical results for ALjGai-^N and In^Gai-^N, 
respectively, and the solid lines are obtained by fitting 
the results with Eq. (Q]). 

After comparing our results with the experimental and 
theoretical works found in literature for both alloys we 
conclude that: 1) Our results for the binary compounds 
are in agreement with the published data, both exper- 
imental and theoretical. 2) For the ALjGai-^N alloy, 
the structural parameters calculated in this work are in 
agreement with those reported in Ref. [131 ] . For this al- 
loy, the bowing parameters for a, c, and u, have positive 
values, indicating a downward bowing. This is in agree- 
ment with experimental and theoretical results reported 
previously. 3) For the In^Gai-^N alloy, the bowing pa- 
rameter of the a and c lattice constants have a nega- 
tive value, while the internal parameter, u, has a positive 
value. Reported theoretical calculations do not give any 
information about the bowing for the u parameter and it 
is not possible to obtain it from experimental measure- 
ments. To the best of our knowledge, this is the first time 
that the value of this parameter has been reported. 



B. Electronic Structure for the ALGai and 
IrixGai-jN alloys 

Before presenting our results of the electronic struc- 
ture, we provide a summary of some representative the- 
oretical and experimental results reported in literature. 

1. Al x Ga l ^ x N 

Hagan et al. 33] and Baranov et al. [34[ were the 
first to demonstrate experimentally the existence of the 
ALGai_2,N alloy. Many other groups have measured its 
lattice constant, c, the optical bowing parameter, <5, and 
the energy gap as a function of the concentration. The 
magnitude of the optical bowing parameter accounts for 
the deviation of the band gap from the linear depen- 
dence. Using MBE, S. Yoshida et al. [l3| measured the 
lattice constant and the band gap for the entire interval 
of concentrations (0 < x < 1). Comparing our results for 
the lattice constant with those reported by this author, 
there is a very good agreement over the entire interval of 
concentrations. Although our calculated band gap value 
shows the experimental trend reported by Yoshida et al. 
[HI], we obtain slightly different values as can be seen 
in Fig. [3] However, other reported values for the band 
gap are well reproduced in our calculations (see Fig. [3]). 
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There are many experimental reports for this system ob- 
tained by different techniques and for different values of 
the concentration. In all these references there is not a 
general agreement concerning if the fundamental param- 
eters, i.e. the lattice constants and the band gap value, 
follow Vegard's law. Positive, negative, or small values 
of the optical bowing parameter can be found through- 
out the literature 0, [H, Q3, G3, [H, |36|. Theoretical 
results have been reported using the kp method, the 
semi-empirical pseudopotentials method, ab initio LDA, 
DFT-LDA using molecular dynamics, and plane waves 
pseudopotentials using DFT-LDA. 

2. In x Ga X - x N 

The first In^Gai-^N alloy with a high degree of or- 
dering in layers grown on sapphire (0001) using metal- 
organic chemical vapor deposition (MOCVD) were ob- 
tained by Ruterana et al. [45| ■ Samples grown with dif- 
ferent concentrations and using different techniques have 
been studied and authors report various values for the 
optical bowing parameter. In Table IIIII we summarize 
the theoretical and experimental results obtained from 
the literature for both alloys along with our calculated 
values obtained for the optical bowing parameter. It can 
be seen from the table that there is a large discrepancy 
between the reported experimental and theoretical val- 
ues. For theoretical calculations, most of the reported 
values are greater than one and they show a larger scat- 
tering of the numerical values, especially for the case of 
the In^Gai-zN alloy. 

In Table ITVl we show the obtained results for the energy 
gap of the ALjGai-zN and In x Gai_ x N alloys from our 
ab initio calculations for different values of the concentra- 
tion, x. These values correspond to 0, 25, 50, 75 and 100 
% Al and In substitution, respectively. The plots corre- 
sponding to these data are displayed in Fig. [31 It can be 
seen that when the concentration of Al is increased the 
ALjGai-^N alloy shows a nearly linear dependence. On 
the other hand, the In^Gai-^N alloy shows a non-linear 
dependence when we increase the In concentration. For 
both alloys we fit the gaps obtained by our calculations to 
a non-linear dependence using the quadratic phenomeno- 
logical function: 

E g {x) = xE 9! a + (1 - x)Eg tB - Sx(l - x); (2) 

where E g> a and E ffi s corresponds to the gap of the A1N 
(InN) and GaN for the Al x Gai_ x N (In x Gai_ x N) alloy. 
Substituting the values of Table llVl into Eq. ^ we found 
6 = 0.3185 and 0.9990 for Al^Gai-^N and In^Gai^N, 
respectively. The solid lines in Fig. 3 correspond to the 
non-linear fitting given by Eq. @. It can be seen from 
this figure that there is a clear non-linear dependence for 
the In^Gai-^N alloy. 

Concerning the electronic structure of these alloys, af- 
ter comparing our results with those reported in litera- 
ture we conclude that: 1) For both alloys, most of the 



experimental results consider concentrations lower than 
x=0.5. This could account for the scattered values re- 
ported for the bowing parameter; however, experiments 
considering the whole interval of concentrations are bet- 
ter reproduced by our calculations. 2) The obtained value 
in this work of S = 0.3185 for the bowing parameter of 
the ALjGai-^N alloy is in agreement with most of the 
experimental and theoretical results reported previously. 
Although many authors claim that the dependence on 
concentration of the band gap parameter of this alloy 
should follow a Vegard's-type law (S = 0.0), it has be- 
come more accepted that there is a small deviation in 
the linear dependence. 3) For the In^Gai-^N alloy the 
reported results show a large scattering, especially in the 
experimental data. This has been partially explained in 
the literature as due to an inaccurate determination of 
the concentration. The quality of the samples and the 
measurement technique also plays an important role in 
the determination of the optical bowing parameter. For 
this alloy, the scattering in the theoretical results re- 
ported in the literature is lower. Our calculations deter- 
mined a value of 6 = 0.9990, which is in good agreement 
with those reported previously. 

Figure [3] shows a comparison of our theoretical results 
with the reported experimental results for the band gap 
energy as a function of concentration for the Al^Gai-^N 
and In^Gai-^N alloys. It can be seen from the figure 
that there is excellent agreement between our results and 
those reported, especially for the Al x Gai_ x N alloy. Most 
of the reported results for the In^Gai-^N alloy are given 
for low concentrations of the In impurity. Nevertheless, 
our theoretical results follow the tendency demonstrated 
by experimental and theoretical results reported previ- 
ously. 

C. Calculated effective masses 

As a by-product of our electronic band structure cal- 
culations, it is easy to compute the curvature of the min- 
imum of the conduction band as well as the maximum of 
the valence band in the vicinity of the T— point. From 
these values the effective masses of the electrons and holes 
can be obtained. 

At the T— point the s— like conduction band effective 
mass can be obtained through a simple parabolic fit using 
the definition of the effective mass as the second deriva- 
tive of the energy band with respect to the wave vector, 
k, via: 

m _ ^mo_&^E_ 

m* ~ h 2 dk 2 ( ' 

where m* is the conduction electron effective mass and 
m is the free electron mass. 

The valence band states at the T— point are derived 
from p— bonding states and for the wurtzite crystals these 
states are not as symmetric as the conduction band. 
However, we can calculate the curvature of the valence 
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band maximum using the following approach: if the spin- 
orbit interaction were neglected, the top of the valence 
band would have a parabolic behavior. This implies that 
the highest valence bands are parabolic in the vicinity of 
the T— point. In this work, all the systems studied satisfy 
this parabolic condition of the maximum of the valence 
band at the T— point Within this approach, and by 
using the appropriate expression of Eq. ^ (using a plus 
sign instead of the minus sign in the prefactor), we have 
calculated the effective masses of the heavy holes at the 
T— point. 

Tabic IVl shows our calculated effective masses for the 
binary compounds, A1N, GaN, and InN, as well as for 
their related alloys. This table also includes theoretical 
and experimental values that were reported in literature, 
for comparison. 

For the binary compounds, we conclude that our cal- 
culated effective masses are in the range of most of the 
reported theoretical and experimental values. Moreover, 
we obtain that our calculated values are in good agree- 
ment with the values recommended by Vurgaftman and 
Meyer 0. 

We report in the same table our calculated effective 
masses for the alloys studied in this work. The calcu- 
lated value for the conduction electron effective mass for 
Al^Gai-zN increases as a function of the Al concentra- 
tion. It should be noted that although we could not find 
published data for this parameter in the wurtzite phase, 
we did find calculated values for the electron effective 
mass of the alloy in the cubic phase [47[ ■ The behavior 
for our calculations and the reported cubic phase values 
show the same trend. The same behavior was also ob- 
tained for the heavy hole effective mass. 

The calculated effective masses for In^Gai-^N show a 
non monotonic behavior as a function of the In concen- 
tration as seen in Table [V] Comparison with empirical 
pseudopotential calculations for the conduction electron 
effective mass in the cubic phase of this alloy [HI, re- 
ported for intermediate values of the In impurity, shows 
that we obtain similar behavior for our calculated effec- 
tive masses. 



D. Zone center phonon calculation 

Using the frozen phonon approach, we can compute the 
zone center phonon modes, Ai (TO) and Ei (TO), by con- 
sidering the change of the total energy as a function of the 
displacement, u, of the atoms in the unit cell from their 
equilibrium positions. In particular, we can compute the 
Ai (TO) mode by considering small displacements of the 
atoms along the optical axis (the c-axis) of the wurtzite 
phase. Then, as we are at the minimum of the total 
energy of the system, the perturbation around this mini- 
mum allows us to use the next parabolic approach to the 
total energy as a function of the displacement: 

E To tai(u) = E + uE 1 + u 2 E 2 = mu; (4) 



Where Ei (i = 0, 1, 2) are fit parameters, [i is the re- 
duced mass, and lo is the frequency. With this approach, 
our calculated value for the A\(TO) mode for the bi- 
nary compounds shows good agreement with most of the 
experimental reports as shown in Table IVII Figure [4^ 
shows our calculated Ai (TO) mode for the Al x Gai_ x N 
alloy, which is compared with the experimental data of 
Refs. 159 . |60J], as well as the theoretical calculations of 
Ref. |61| . Although, we obtain slightly different values 
over the range of Al concentration, we observe that the 
figure depicts the same trend shown by the experimental 
reports. Figure UJd shows our calculated Ai (TO) mode 
for the In^Gai-^N alloy and compares them with theo- 
retical calculations reported in Ref. [61] . As in the pre- 
vious case, it can be seen that our calculated values for 
the Ai (TO) mode reproduce the trend reported in Ref. 
[6l| . In this way we show that our calculations properly 
reproduce the reported values of the Ai (TO) mode for 
the alloys and their binary parent compounds studied in 
this work. 



IV. CONCLUSIONS 

We have calculated the structural and electronic prop- 
erties of wurtzite A1N, GaN, InN, and their related al- 
loys, ALrGai-zN and In^Gai-^N. We found that, for 
both alloys, their structural parameters as a function of 
the concentration, x, do not follow Vegard's law. We 
observe that for the Al^Gai-^N alloy, the a, c, and u 
parameters have a positive bowing parameter, of which, 
the lattice constant, c, demonstrates the largest value. 
On the other hand, when the concentration increases in 
the In^Gai-zN alloy, the bowing parameter for the a 
and c lattice constants have a negative value and the 
bowing parameter for the internal parameter, it, remains 
positive. From our results of the electronic band struc- 
ture calculations we obtained the band gap energy as a 
function of the concentration, x, and characterized the 
deviation from the linear behavior by calculating the op- 
tical bowing parameter, 5. We obtain a small optical 
bowing parameter 5 — 0.3185 for the Al^Gai-^N alloy, 
which is consistent with results reported previously. The 
deviation from the linear behavior is more drastic for the 
In-rGai-^N alloy, which has 5 = 0.9990 in accordance 
with most of the experimental results found in litera- 
ture. The effective masses of the systems studied in this 
work were calculated and we found that our calculated 
effective masses for the binary compounds, A1N, GaN, 
and InN, are in the range of the reported values in the 
recent literature. To the best of our knowledge, this is 
the first time that the deviation parameter for the in- 
ternal parameter, e u , as well as the effective masses for 
these alloys have been reported. Finally, using the frozen 
phonon approach we have computed the Ai(TO) mode 
for the different systems studied in this work. From our 
calculations we have found good agreement with exper- 
imental values reported for the binary compounds. For 
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the ternary alloys, our calculations reproduce experimen- 
tal values for Al^Ga^^N, as well as theoretical predic- 
tions for In x Gai_a;N. 
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TABLE I: Structural parameters (in A) for the A1N, GaN, and Al^Gai-^N alloys. The lattice parameters a and c are given in 

A. 





Present Work 


Exp. Results 


Other Calc. 


GaN a 


3.2209 


3.1890", 3.1892 6 
S.1880% 3.19 9 


3.1660 c , 3.189 d , 3.2 C 
3.1800\ 3.1986 1 , 3.17 J , 
3.183* 


c 


5.2368 


5.1850 a , 5.185 
5.18561 / , 5.189 9 


5.1540 c , 5.185 d , 5.2200 e 

5. lags' 1 , 5.2262 1 5.15P, 

5.178* 




u.o / ou 


n "i7P,R a n ^77 9 

U.O / UO , U.OI 1 


n ^77n c n ^7fis <i n i7F,n e 

U.OI IU , U.OI UO , U.O 1 uu 

n ^7(\n h n ^779* n ^p,^ 

U.O(UU , U.O 1 , U.O 1 UO 


A In ncflan tcIM A 
- rii U.2o ^"^U. 75 1 ' <* 


3 2059 


Ref m rn 


3.163* 




5.1338 




5.137* 




n ^781 

U.OI 01 






Alo.50Gao.50N a 


3.1719 


See Ref. [H, 


3.139* 


c 


5.1012 




5.085* 


u 


0.3790 






Alo.75Gao.25N a 


3.1601 




3.098* 


c 


5.0870 


See Ref. [HGJ] 


4.990* 


u 


0.3791 






A1N a 


3.1411 


3.1120 a , 3. II 9 
3.1106' 


3.0920 c , 3.084 d , 3.1 e 
3.0610\ 3.10954% 3.098 J , 
3.076* 


c 


5.0268 


4.9820 a , 4.98 9 
4.9795' 


4.954 c , 4.9948 d , 5.010 e 
4.8976\ 4.9939\ 4.9599 J ', 
4.935* 


u 


0.3805 


0.3819 a , 0.3821 9 


0.3821 c , 0.3814 d , 0.3800 e 
0.3820\ 0.3819% 0.3819- 7 ' 



a Ref. 
'Ref. 
c Ref. 
d Ref. 
e Ref. 
/Ref. 
9 Ref. 
h Ref. 
'Ref. 
^Ref. 
*Ref. 
; Ref. 
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TABLE II: Structural parameters 

A. 



(in A) for the GaN, InN and In^Gai-^N alloys. The lattice parameters a and c are given in 



"Ref. 
6 Ref. 
c Ref. 
d Ref. 
e Ref. 
/Ref. 
9Ref. 
fe Ref. 
*Ref. 
JRef. 
*Ref. W 
! Ref. 
™Ref. 
"Rcf. 



30 
31 
32 







JT I fcJDfcJH L VVUIK 


JZjXJJ . rVcolllLo 


W Lllfcil KJOjIKj. 


\JCLL\ 


M 


Q 990Q 


o. ioyu , o. ioyz 


o. iuuu , o. ioy , o.zuuu 








O.lOOU - O.lc/ 


3 1 R00 h ^ 1 QSfi* 3 1 7 J 
o.iouu , o.iyou - o.i/ 




C 


^ 9^8 


O.loOU , O.loO 


^ 1 ^AC\ C ^ 1 ^ 990n e 
O.lO^U , O.loO , O.ZzUU 








^ 1 8^81 ^ ^ 1 RQ g 

O . ±00\J -L - u.lOc/ 


u-ioyo ; o.ZiZiUZi . j.i ji 






U.O ( OU 


^7(SS a 177 9 

U.O / Uu , U.O 1 I 


^77n c ^7ftR d *37fifl e 

U.O( (U , U.OIUO j U.O I uu 










o ^7fin^ n ^779^ n ^7fi&? 

U.O/UU , U.O I I — j U.O I uo 


Tffn ncHfln "re 1ST 




t J • iJj£j£7CJ 










5.3987 










fl ^7Q1 

U.O ( U± 






1110.50 VI *±0.50 




3.4128 










5.5257 










n ^7Q9 
u.o ( 






1110. 75 v* a-0.25^ 












c 


O.DOoo 










u.o ( yu 






InN 


a 


3.5440 


3.5365*, 3.5378' 


3.520 c , 3.501% 3.480 e 








3.548 m , 3.540" 


3.525% 3.614% 3.546 J 




c 


5.7228 


5.7039% 5.7033' 


5.675 c , 5.669% 5.64 e 








5.76 m , 5.705™ 


5.68583% 5.8836% 5.7162 J 




u 


0.3806 




0.3799% 0.3784% 0.378 e 










0.379% 0.37929% 0.379^ 
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TABLE III: Experimental and theoretical results for the optical bowing parameter 8 for ALjGai-ajN and In x Gai_ x N. 



Al^Gai-sN 



Experimental 


S [eV] 


Theoretical 


8 [eV] 


@] Y. Koide, for < x < 0.4 using MOVPE 


w 1 


[ll| M. Goano using pseudopotentials and VCA 


0.069 


[ll] S. Yoshida, < x < 1 using MBE at 700° 


w 


Z. Dridi using FP-LAPW LDA and VCA 


0.710 


[3 F. Yun, for < x < 1 using MBE at 600-670° 


1 


[13 D. K. Wickenden,0 < x < 0.4 using MOCVD 


fa 


[H M. A. Khan, < x < 1 using MOCVD at 


w 


[37J K. Chen using ab initio Molecular dynamics 


1.40 


915° 








[H T. Takeuchi, for < x < 0.25 using MOVPE 


w 1 


[38(] S. K. Pugh using first principles k p method 


fa 


[39l] 0. Ambacher Review article 


1 


Value calculated in this work 


0.3185 


In^Gai-^N 


Experimental 


5 [eV] 


Theoretical 


8 [eV] 


[H V. Yu. Davydov et al. for 0.36 < x < 1 


2.50 


[TH] M. Goano et al. 


1.05 


HI] T. Takeuchi et al. 


3.20 


QJ| Z. Dridi al. 


1.70 


[U S. Nakamura et al, for < x < 0.4 


1.00 


[43| C. Caetano oi., using FP and LDA 


1.44 


gl[ M. D. McCluskey et al, for < x < 0.12 


3.5 


[11 A. F. Wright et al. 


0.1677 


[11] C. Wetzel et al, for < x > 0.2 using 


2.6 


Value calculated in this work 


0.999 


MOVPE 








[3 J. Wu et al, for < x < 0.5 using MBE 


1.4 







TABLE IV: Band gap energy (eV) for the Al x Gai_ x N and In x Gai_ x N alloys as a function of the concentration, x, obtained in 
the present work. 



System 


E 9ap (eV) 


GaN 


1.768 


Alo.25Gao.75N 


2.319 


Alo.50Gao.50N 


2.830 


Alo.75Gao.25N 


3.5123 


A1N 


4.027 


GaN 


1.768 


Ino.25Gao.75N 


1.088 


Ino.50Gao.50N 


0.871 


Ino.75Gao.25N 


0.489 


InN 


0.299 
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TABLE V: The electron and hole effective masses for A1N, GaN, InN and their alloys, m* denote the average masses. The 
average effective mass can be obtained using m* = [mp_, M m,p^ K mjl_ (A ] 1 ' /3 . Where, m x and m" denote the k direccion 
dependent masses perpendicular and parallel to the c axis, respectively. All values are in units of a free-electron mass mo. For 
a recommended set of values for the binary compounds see Ref. |U. 



a Ref. 
'Ref. 
c Ref. 
d Ref. 
e Ref. 
•faef. 
9Ref. 
h Ref. 
^ef. 
•»Ref. 
fc Ref. 
'Ref. 
m Ref. 
"Ref. 



49 
58 
29 
56 
57 
50 
51 
53 
52 
54 
57] 
62] 
63] 
64] 







Present work Other results 






m e 


TUh 


m e mh 


A1N 


m x 
mil 
m* 


0.3012 
0.2847 
0.2956 


4.3243 
0.2427 
1.6528 


0.30 a 0.33 6 0.33 c 
0.32° 0.32 6 0.32 c 
0.31 d 0.48 c 


4.35 c 
0.28 c 


Alo.75Gao.25N 


m x 
m" 
m* 


0.2682 
0.2913 
0.2749 


3.5247 
0.9325 
1.4743 






Alo.50Gao.50N 


m x 
m" 
m* 


0.2330 
0.2411 
0.2427 


2.3039 
0.5345 
1.2961 






Alo.25Gao.75N 


m x 
m" 

* 

m 


0.2020 
0.1958 
0.2000 


3.1232 
0.1741 
1.1969 






GaN 


m x 
m" 
m* 


0.1491 
0.1803 
0.1692 


2.1072 
2.1048 
2.1412 


0.20 a 0.21 6 0.22 c 
0.20° 0.19 6 0.20 c 
0.22 7 0.23 9 0.20 h 
0.18 J 0.20* 0.24 ! 
0.215" 1 


0.39 c 
2.04 c 
2.2* 


Ino.25Gao.75N 


m x 
m" 
m* 


0.1069 
0.0968 
0.1035 


2.5409 
2.4447 
2.4953 






Ino.50Gao.50N 


m x 
m" 
m* 


0.1182 
0.0921 
0.1025 


2.7958 
2.6116 
2.8371 






Ino.75Gao.25N 


m x 
m" 
m* 


0.0717 
0.0623 
0.0781 


2.0304 
2.1486 
2.2179 






InN 


m x 
m" 
m* 


0.1299 
0.0892 
0.1146 


1.9096 
2.2051 
2.0301 


0.07 a 0.068 6 0.07 c 
0.068" 

0.07 a 0.065 6 0.06 c 
0.072" 


2.967" 
2.566" 
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TABLE VI: Calculated Ai(TO) mode for the different sys- 
tems studied in this work. Comparison with other calcula- 
tions and experimental data (all values in cm -1 ). 



System 


Present work 


Other results Experimental 


GaN 


516 


537° 545 6 


533.8 C 533. 5 d 


Al 25 Ga75N 


529 






Al 50 Ga 5 oN 


538 






Al 2 5Ga 25 N 


548 






A1N 


565 


649 a 615 c 619 e 


613.8 C 


In 2 5Ga 7 5N 


507 






In 5 oGa 50 N 


491 






In 7 5Ga 25 N 


475 






InN 


466 


443 ; 


447 9 440 / 



a Ref. 
^Ref. 
c Ref. 
d Ref. 
e Ref. 
/Ref. 
»Ref. 



65] 
66] 
67] 
68] 
69] 
70] 
711 
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0.00 0.25 0.50 0.75 1.000.00 0.25 0.50 0.75 1.00 

Al concentration x In concentration x 

FIG. 1: (a) Lattice constants a and c for the ALGai-^N alloy as a function of the Aluminium composition, x. Open circles 
show the experimental results from Ref (b) Lattice constants a and c for the In^Gai-^N alloy as a function of the Indium 
composition, x. 
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0.00 0.25 0.50 0.75 1.00 



Concentration x 



FIG. 2: Internal parameter u for the Al^Gai-^N and In^Gai-^N alloys a function of the Aluminium and Indium composition, 
x, respectively. 
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0.0 1 1 1 1 1 1 1 1 1 

0.00 0.25 0.50 0.75 1.00 

Concentration x 



FIG. 3: Variation of the band gap for the Al^Gai-^N and In^Gai-^N alloys as a function of the concentration, x. We present 
the experimental results from Table 3 and the solid lines represent the proposed adjustment to our results (solid points) 72]. 
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FIG. 4: Ai (TO) phonon dependence with Al concentration for Al^Gai-^N alloy (upper figure), and as a function of the In 
concentration for In^Gai-^N (lower figure). We present our calculations (full circles), the experimental results of Ref. [HtJ 
(open circles with plus), and Ref. [6(J (open circles). Solid line depicts the calculated values from Ref. [6ll |. 



